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The conducting and thermodynamic properties of ballistic metallic nanocontacts with smooth 
shapes are investigated. All properties are related to the electronic scattering matrix, which 
is evaluated in the WKB approximation for independent electrons and in the self-consistent 
Hartree approximation for interacting electrons. Mesoscopic oscillations of order InN in the 
cohesive force and of order e in the contact charge are predicted when a metallic nanocontact 
is pulled apart, which are synchronized with quantized jumps in the conductance. 



1 Introduction 

In a seminal experiment published in 1996, Rubio, Agra'ft, and Vieira 1 measured simultaneously 
the electrical conductance and cohesive force of Au-Au contacts with diameters ranging from 
several Angstroms to several nanometers. They observed oscillations in the cohesive force of 
order InN as the contact was pulled apart, which were synchronized with jumps of order Go = 
2e 2 jh in the conductance G. Similar results were obtained independently by Stalder and Diirig? 
In a previous article, 3 it was argued that the mesoscopic force oscillations, like the corresponding 
conductance steps, can be understood by considering the nanocontact as a waveguide for the 
conduction electrons (which are responsible for both conduction and cohesion in simple metals): 
each quantized mode transmitted through the contact contributes 2e 2 /h to the conductance, 4 
and a force of order Ep/Xp to the cohesion, 3 where Xp is the de Broglie wavelength at the 
Fermi energy Ep. Recent experiments 5 suggest that mode quantization plays an important 
role in the cohesion of Na nanocontacts for G < 120Go- In Ref. 3 and several subsequent 
works, 6-9 an independent-electron model was used to describe the cohesion and conductance 
of metallic nanocontacts within the grand canonical ensemble. Recently, several attempts have 
been made to include the effects of electron-electron interactions in the extended Thomas- Fermi 
approximation 10 and in the local density approximation; 11 however, these calculations have 
utilized the canonical ensemble, which is not appropriate for an open mesoscopic system, such 
as a metallic nanocontact. 

In this article, we show how to treat electron-electron interactions self-consistently in open 
mesoscopic systems via a generalization of the scattering-matrix approach of Ref. 3. In Sec. 
2, we describe the scattering-matrix formalism for independent electrons, and show that the 
bare mesoscopic charge oscillations induced by the quantum confinement are small, of order e in 
nanometer scale contacts. In Sec. 3, we derive the self-consistent Hartree approximation using 
the electronic scattering matrix, and obtain integral equations for the screened charge and grand 
canonical potential within linear response. The Hartree equations are solved for the case of a 
discrete-potential model, and it is shown that the Hartree correction to the cohesive force is 
small, justifying the use of the independent-electron model. 3 > 6-9 
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2 Independent electron model of a metallic nanocontact 



We first model the metallic nanocontact as a system of noninteracting electrons confined along 
the z-axis by a potential E/"o(x) with a saddle-point at x = 0, and assume Uq varies slowly, 
i.e., dUo/dz <C Ep/Xp Vx = (x,y,z). Neglecting small terms of order dUo/dz, the Schrodinger 
equation is separable into a series of one-dimensional scattering problems with effective potential 
barriers e u (z), the eigenvalues of the two-dimensional Schrodinger equation in planes of constant 
z. The scattering matrix is then diagonal, and for a system with inversion symmetry about z = 0, 
it is found in an extended WKB approximation to be 3 



v(fi) = i„^| T \ /2 {E) ), (1) 



where 



e„(E) = (2m/h 2 ) 1 / 2 f dz[E- Ev(z)] 1 ' 2 (2) 

Js v (z)<E 

is the scattering phase shift, and R V {E) and T V {E) are the reflection and transmission proba- 
bilities for mode v, which satisfy R U (E) + T U (E) = 1. The transmission probability for mode v 
may be expressed 12 as T^(E) = 1 + ^(E 1 ) -1 , where 



ln%(E) 



-2{2m/% 2 ) 1 ' 2 f dz [e v {z) - E} 1 / 2 , E < e„(0), 

j£^(z)>£; 

(3) 

2(2m/h 2 ) 1 / 2 f dz[e u (z) + E-2e u (0)} 1 / 2 , E > e u (0). 

Je v (z)>2e v {0)-E 



Given the scattering matrix (l)-(3), the electrical conductance of the nanocontact may be 
obtained from the Landauer formula, 4 and all equilibrium thermodynamic properties of the 
system may be obtained from the density of states g(E), 



(4) 



dE ' 

The grand canonical potential, which governs the energetics of the nanocontact is 3 

n = dEg(E){E -E F ) = - 8 ^j:f dz (l - ^l) 3 ' 2 . (5) 

JO oAf ~ Je v {z)<E F \ &F / 

The total number of electrons in the nanocontact is 

N-= [ EF dE 9 (E) = ±^[ dz (l - £ -^) 1/2 . (6) 

JO *F „ Js v (z)<E F \ E F J 

Eqs. (4), (5), and (6) are given at zero temperature for simplicity. When the nanocontact is 
elongated, the contact area decreases, steepening the saddle in the confining potential, and hence 
increasing the mode energies e v (z). The cohesive force is given by F = —80,/dL. 



2.1 Hard-wall confinement 

We now specialize to the case where electrons are confined within the nanocontact by impene- 
trable walls, with Uq = on the interior. In this case, the condition for the validity of Eqs. (5) 
and (6) is \dD(z) / dz\ <C 1, where D(z) is the diameter of the contact. 3 In order to separate the 
macroscopic contributions to Q and iV_ from the quantum fluctuations, it is useful 6 to express 



the sum over transverse modes v in Eqs. (5) and (6) in terms of an integral over a local transverse 
density of states ^2 U J dz (• • •) — > J dE J dz dNj_(E, z)/dE (• • •). The total number of transverse 
modes N±(E, z) with e v (z) < E may be determined with the aid of the Weyl expansion 6 

N^E,,) = ^-H^ + \ + tN x (E,,), (7) 

where A(z) and dA(z) are the cross-sectional area, and circumference, respectively, of the contact 
as function of z, and 5N±(E,z) is a fluctuating correction determined by the spectrum of 
quantized transverse modes, whose average is zero. Using Eq. (7), Eqs. (5) and (6) may be 
rewritten as 

n 2klV klS 2k F L 5n A7 k%V klS k F L rA7 

— = ^rr + — 1 , N- = -V £— + — h 6N-, (8) 
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where V is the volume of the nanocontact, S is its surface area, L is its length, and 5Q and 5iV_ 
are oscillatory quantum corrections whose average is zero. 

During elongation, it is assumed that the positive background charge density remains con- 
stant, i.e., that the system is incompressible. Hard- wall boundary conditions impose a node of 
the electron wavefunction at the boundary, and lead to a depletion of electrons within a distance 
of order Xf from the boundary [second term on the right of Eq. (8)]. To reconcile the constraint 
of incompressibility and the Dirichlet boundary conditions, we assume that the shape D{z) of 
the nanocontact evolves during elongation in such a way that the macroscopic part iV_ of the 
electron number of the system [first three terms on the right of Eq. (8)] remains constant, 13 and 
is neutralized by an equal and opposite positive background charge. The net charge imbalance 
on the nanocontact (neglecting screening) is thus 5Q$ = —e5N_, which we show to be small. 
Differentiating Q with respect to L with the constraint iV_ = const., and using Eq. (8), one 
finds 



F ~~dL 



where <tv = Epkp /16ir is the surface energy 3 of a noninteracting electron gas at fixed V and 
AFtop = AEp/9Xp is the mean mesoscopic suppression 6 of the cohesive force at fixed V. The 
reduction of the surface energy by a factor of 5 has been discussed by Lang. 1 Importantly, 
since the constraint iV_ = const, differs from the constraint V = const, only by terms of order 
(fcF-Dmin) -1 , the mesoscopic fluctuations 5F and 5N- are quite insensitive to the choice of 
constraint. 

Figure 1 shows the electrical conductance G of a metallic nanocontact, calculated from 
Eqs. (3) and (4), and its cohesive force F, calculated from Eqs. (5) and (9). The conductance 
decreases in steps of 2e 2 /h and 4e 2 /h (depending on whether the transverse mode involved is 
singly or doubly degenerate) as the contact is elongated, and the force exhibits oscillations of 
amplitude ~ Ep/Xp, which are synchronized with the conductance steps. The overall magnitude 
of the force is reduced relative to that obtained in a previous calculation 3 using a constraint of 
constant volume, but the force oscillations 5F are virtually identical to those obtained in Ref. 3. 
Noting that Ep/Xp ~ 0.76nN for Na, one finds that the force shown in Fig. 1(b) is also virtually 
identical to that calculated for a Na nanowire by Yannouleas et al., 10 suggesting that the more 
realistic jellium confinement potential they utilized does not significantly modify the cohesion. 

The net charge 5Qq on the nanowire (neglecting screening) is also shown in Fig. 1. Interest- 
ingly, we find that 5Qo/Lq ~ e5F/Ep, where Lq is the initial length of the (initially cylindrical) 
nanocontact. This result is consistent with the result of Kassubek et al. for an abrupt geometry, 
and will be discussed in more detail elsewhere. Here it suffices to note that since \5Qq\ <C eiVL, 
the electron-electron interactions in the nanowire may be treated using linear-response theory. 
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Figure 1: Properties of an axially symmetric nanocontact with diameter D(z) = D m in + (£>max — D m m)(2z/Z/) 2 as a 
function of the relative elongation AL/Lq, where fcpD maj( = 24 and ftp-io = 40: (a) Electrical conductance G; (b) 
cohesive force F in the independent-electron model (solid curve) and in the self-consistent Hartree approximation 
(points); (c) mesoscopic charge oscillations 8Qo (neglecting screening) and force oscillations SF; (d) self-consistent 
Hartree correction to the cohesive force in the constant interaction model [Eq. (15)]. 



3 Interacting electron model 

3.1 Self- consistent Hartree approximation 

In an interacting model, the confinement potential of the nanocontact is determined self-consistently 
via the Coulomb interaction: Uq{x) — > U(x) = J a?y V(x — y)6p(y) (Hartree approximation), 
where <5p(x) = P-(x) + p+(x) is the net charge density, V(x) is the Coulomb interaction, and 
p+(x) describes the fixed positive background. The grand canonical potential becomes 

n = n [u(*)] -\j<?x [ P _(x) - P+ (x)] i/(x), (io) 

where f2o[^(x)] is given by Eq. (5) with t/o( x ) replaced by f/(x). The second term in Eq. (10) 
corrects for the double counting of electron-electron interactions in f2o[J7(x)]. The charge density 
of the conduction electrons is 

3.2 Linear response 

Motivated by the smallness of the bare charge fluctuations 5Qo = J d 3 x5po(x.), we expand Eq. 
(11) to linear order in f/(x), obtaining a self-consistent equation for the screened charge density 



/ 



d 3 y e(x, y)5p(y) = <5p (x), e(x, y) = S(x - y) - J d 3 z 



<Jp(x) 



5U(z) 



V(z-y), (12) 
u=o 



where e(x, y) is the nonlocal Hartree dielectric function and 

<5p(x) 1 [ Ef \ <fS f 6S ot 5 2 S 



<Stf(y) 2W0 I «%) <5C/(x) + 6 ^(y)«JC/(x) " H ' C ' f ' (13) 



Expanding Eq. (10) in powers of U(x) and using Eq. (12), one finds 

n = n [0] + ±Jd 3 xJ d 3 y 6p (x.)V(x,y)6p (y) + O(5p 3 ), (14) 

where the screened interaction is given by V(x, y) = / d 3 z V(x — z)e _1 (z,y) and the inverse 
dielectric function e" 1 satisfies / d 3 z e _1 (x, z)e(z, y) = / d 3 z e(x, z)e _1 (z, y) = <5(x — y). Eq. 
(14) is a familiar result of linear-response theory; the new element here is that the dielectric 
function is calculated from the scattering matrix [Eqs. (12) and (13)] and is nonlocal due to 
quantum confinement. 

3.3 Constant interaction model 

A simple case for which Eqs. (12) and (14) can be solved easily is the so-called constant inter- 
action model, V(x, y) = C _1 if x, y G V and otherwise. C plays the role of the capacitance of 
the nanocontact to its environment, and Eqs. (12) and (14) become: 

The result (15) confirms a previous conjecture 7 based on simple physical arguments. An upper 
bound on the magnitude of AQh = 0, — Qo[0] is obtained by setting C = 0, which enforces 
global charge neutrality. The Hartree correction to the cohesive force AFh = —dAQn/dL is 
plotted in Fig. 1 for C = 0. AFh is over an order of magnitude smaller than 5F, and gives 
an indiscernible correction to the total force, suggesting that the independent electron model of 
metallic nanocohesion 3 ' 6-9 is well justified. 
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